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NON-VARIATIONAL APPROXIMATION OF DISCRETE 
EIGENVALUES OF SELF-ADJOINT OPERATORS 

LYONELL BOULTON^ 


Abstract. We establish sufficiency conditions in order to achieve 
approximation to discrete eigenvalues of self-adjoint operators in 
the second-order projection method suggested recently by Levitin 
and Shargorodsky, CSl We find explicit estimates for the eigen¬ 
value error and study in detail two concrete model examples. Our 
results show that, unlike the majority of the standard methods, 
second-order projection strategies combine non-pollution and ap¬ 
proximation at a very high level of generality. 


1. Introduction 

Let M be a self-adjoint operator acting on a dense domain, DomM, 
in a separable Hilbert space Ti. Let the spectrum of M be 
split into isolated eigenvalues of finite multiplicity (discrete spectrum) 
and degenerate points (essential spectrum), in symbols 
Spec M = SpecjiscAf U SpeCgsgM. In order to approximate discrete 
eigenvalues A of M, we may consider using the following projection 
method: choose an orthonormal basis of 7i, C DomM, and 

hnd the spectrum of large matrices M„ resulting from compressing M 
to the hnite dimensional subspaces £„ = spanj^i, 

The successful outcomes of this strategy is illustrated by the well 
known Rayleigh-Ritz theorem for the approximation of variational eigen¬ 
values, see e.g. m Theorem XIIL4]. Assume that M has a non¬ 
degenerate ground eigenvalue A = min [Spec M] < min[SpeCess Af]. Let 
n„ be the orthogonal projection onto £„ and M„ := n„M|£„. If 

(1) lim ||M„n„'0 — A'011 = 0 for all = Xip, 

n—>co 

then the first eigenvalue of Mn converges from above to A (i.e. ap¬ 
proximation is achieved) and the second eigenvalue of M„, counting 
multiplicity, can not be smaller than min[SpecM\ {A}] (i.e. no chance 
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of spurious eigenvalues). Conditions such as ((H) are useful in applica¬ 
tions because they can sometimes be verihed on abstract grounds, CHI 
Ex. XIII.2.1 and XIII.2.2]. Similar results hold for any other variational 
eigenvalue (i.e. those below the minimum of the essential spectrum or 
above the maximum of this set) simple or otherwise. 

Although this technique is proven to be valuable in the study of 
spectral properties of self-adjoint operators, think e.g. of Galerkin ap¬ 
proximation, two main difficulties are in place when we extend it to 
the analysis of non-variational eigenvalues: 

- Lack of approximation: no eigenvalue of is close to A in the 
limit n —>■ cx), and 

- Spectral pollution: there are eigenvalues Zn of which appear 
to converge to some p not in the spectrum of M. 

In fact any /i G conv[SpeCeggM] \ SpecM is potentially an accumula¬ 
tion point of SpecM„ in the limit n oo. This is known to occur 
in several important applications such as elasticity, electromagnetism 
and hydrodynamics, see e.g. ID, M, Hi and the references therein. 
Widely available commercial packages such as FEMLAB are known to 
produce spectacularly incorrect results even in cases of the simplest 
one-dimensional Stokes’ type systems, ng. Actually one may con¬ 
struct simple examples of bounded M and natural orthonormal basis 
[fk] with the property that every p G conv^peCgg^M] \ SpecM is an 
accumulation point of SpecM„ in the limit n —> oo, see dH and m. 

The standard approach to deal with pollution around non-variational 
eigenvalues, aims at choosing £„ sufficiently close to the spectral sub¬ 
spaces of M in order to capture the eigenfunctions associated to A, see 
e.g. d], 12] and dni- Unfortunately no universal device is known for 
constructing such subspaces. 

A different approach hinges on the spectral theorem and it exploits 
the fact that for C, lying in a gap of SpeCess M, the eigenvalues of C — M 
inside the corresponding shifted gap enclosing the origin, become vari¬ 
ational eigenvalues of {( — M)^. In dHj Levitin and Shargorodsky 
consider this idea in a concrete manner. They propose complement¬ 
ing the usual projection method with a posteriori estimates found by 
computing the eigenvalues of the matrix polynomial 

(2) Pn{z) := - My\Cn = z^- 2MnZ + [MX. 

Usually Spec Pn contains non-real points (see the dehnitions in Sec¬ 
tion 2). Their result is based on a remarkable property |151 Theo¬ 
rem 2.6] guaranteing that for 2 ; G SpecP„, 

(3) [Rez—|Imz|, Rez-|-|Imz|] fl SpecM 7 ^ 0 . 
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Hence the strategy would be finding points in SpecPn- If these points 
are close to R, they must also be close to the spectrum of M. 

Numerical evidence found in j2] and ra suggests that approximation 
occurs for concrete operators and natural basis {4>k}- The aim of the 
present paper is to discuss in detail this complementary problem of 
approximation in this “second-order” projection method. Our main 
contribution summarizes in the following statement: i/A G SpeCdiscTf 
and 

(4) lim — A^-^ll = 0, A; = 0,1, 2, whenever Mip = 

n—^■oo 


then there exists Zn G Spec Pn such that Zn ^ X as n ^ oo. See 
Theorems HI and 121 below and also jl]. Hence, discrete eigenvalues sat¬ 
isfying © will always be approached no matter their location relative 
to SpeCessTf. Notice that condition 0 here plays the role of o in the 
“linear” projection method. 


In Section 2 we establish the main theoretical contribution. Theo¬ 
rem^ finds explicit estimates for the eigenvalue error — A| in terms 
of bounds for the difference in the left side of ©• This result depends 
on knowing exactly the entries of the matrices Mn and hence 

its scope is mainly theoretical. Theorem |21 on the other hand, demon¬ 
strates that introducing noise to the entries of Pn{z) within a certain 
tolerance, does not affect approximation. This result is better suited 
for applications. We give precise bounds for this tolerance in terms of 
invariants of the problem. 

Section 3 illustrates the scope of the theoretical results by means 
of two model examples. The first one corresponds to finite rank per¬ 
turbations of multiplication operators by a bounded symbol. We find 
explicit estimates for the error \zn — A| when {cpk} is of Fourier type 
and report on numerical outputs. The second model is unbounded, the 
Schrodinger operator with band gap essential spectrum. We choose 
{cpk} to be the Hermite functions, then find explicit upper bounds for 
eigenvalue error in terms of smoothness properties of the potential. 

In order to make the paper more readable, the proofs of Theorems ^ 
and 121 are to be found in Section 5 while Section 4 introduces and 
describes a convenient notation in order to simplify the presentation 
of most arguments. Being of obvious interest, in Section 6 we discuss 
the question of whether approximation also occurs for the essential 
spectrum. 
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2. Speed of convergence and stability: general results 

We begin this section by fixing some notation. Let A G Specdisc^- 
Throughout this paper {n„} denotes a family of orthogonal projec¬ 
tions onto Ti satisfying dl such that Cn '■= Rann„ C DomM^. The 
sequence 5{n) >0 shall always denote an upper bound for the left side 
of 0) such that 5{n) —>■ 0 as n —>■ oo\ for all "0 G DomM satisfying 
M'lp = X'lp and H-^H = 1, 

(5) - A>|| < (5(n) for n G N and fc = 0,1, 2. 

Below we shall always assume that "0 is an eigenvector associated to A 
normalized by H-^H = 1 . 

Let P{z) = where zisa. complex variable, G and 

det Am 7 ^ 0. The spectrum of P is, by definition, the set of eigenvalues 

SpecP := { 2 ; G C : det P( 2 ;) = 0}. 

The hypothesis det Am 7 ^ 0 ensures that Spec P comprises no more 
than mn finite points. 

We are concerned with the spectrum of the quadratic matrix polyno¬ 
mial Pn{z) given by ©• The determinant of Pn{z) is a scalar polyno¬ 
mial in 2 ; of degree (2 dim £„), so that Spec P„ is a finite set comprising 
at most (2dim£„) different complex numbers. In general these points 
do not intersect the real line, unless £„ contains an eigenfunction of M. 
Furthermore, since det Pn{z) = det Pn(T), M is an axis of symmetry for 
SpecP„. In Section 4 we characterize SpecP„ as the poles of a certain 
positive-valued subharmonic function. Other descriptions of the spec¬ 
trum of matrix polynomials better suited for computation, e.g. as the 
eigenvalues of block matrices, may be found in the literature, see ma 
or nn Part II]. 

The goal of the procedure suggested by Levitin and Shargorodsky in 
na is finding the points in Spec P„ which are close to R. Statement Q 
ensures that these points will necessarily be close to SpecM. A com¬ 
plementary assertion guarantees that discrete eigenvalues satisfying (jlj) 
will always be approached no matter their location relative to SpeCessdL. 
Our first theorem in this direction establishes that the rate of conver¬ 
gence | 2 ;„ — A| —> 0 for G SpecP„ is at least a power 1/2 the rate at 
which the Cn approximate k = 0,1, 2. It is unclear whether the 

claimed power 1/2 is sharp for general choices of Cn- In Section l!0 
below, we provide specific numerical evidence suggesting that for par¬ 
ticular cases the actual eigenvalue error is — A| = 0(5(n)") for 
a ^ 1. 
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Theorem 1. Let A G Spec discAf and assume that © holds. There 
exists 6 > 0 independent of n and Zn G Spec Pn, such that 

( 6 ) \zn — A| < b[S{n)Y^'^, n G N. 

The scope of Theorem ^ is mainly theoretical. Computing Spec Pn 
usually requires estimating the coefficients of Pn{z). Since Pn{z) is 
Hermitian for all z G M, a well known result in the theory of matrix 
polynomials (cf. [TTl Theorem II.2. 6 ]), ensures that we can always hnd 
a factorization 

(7) Pr,{z) = iz-A^){z-A*J, zee, 

where An are square matrices of size dim£„ constant in 2 ;. Clearly 
Spec Pn = (Spec An) U (Spec An). A small perturbation in the coeffi¬ 
cients of Pn will change the coefficients of An- Typically the eigenvalues 
of An are not semi-simple. In fact their condition number might, in 
principle, be large as n is large, forcing Spec Pn to be sensitive to small 
changes in the entries of Mn and [M^]n, see [12]. Our next result 
establishes that this sensitivity can be controlled uniformly in neigh¬ 
bourhoods of SpeCdiscAf. To be more specihe, approximation to a small 
5-neighbourhood of A is always achieved, if we estimate the coefficients 
of Pn within an error smaller than some prescribed tolerance, Wk£. 

Below and elsewhere we shall write p := dist [A, SpecM \ {A}] > 0. 
The norm || ■ || for matrices shall always refer to the uniform operator 
norm. Let £ > 0 and w = (wo,Wi) where Wk > 0. We denote by P^^n, 
the set of sequences of linear matrix polynomials (Qn)^=i such that 
Qn(z) = FnZ — Gn, where Fn and Gn are square matrices constant in 
z of size (dim£„) satisfying ||Gn|| < wqE and ||TA|| < wis. 

Theorem 2. Let A G SpeCdiscAf. Assume that holds and that 
Wo,Wi >0 do not vanish simultaneously. Lett) < 5 < p/4 be fixed. Let 

^ ^ 2(252 + 3p^)[wo + wfip/i + |A|)] ■ 

There exists N > 0 only dependant upon 5, wq and wi, such that for 
all {Qn)ffL, G 

Spec {Pn + Qn) n {5 < — A| < p/4} = 0 and 

Spec {Pn + Qn) n {I 2 ; — A| < 5} 7 ^ 0 for n > W 

Moreover, if we count multiplicities, the number of eigenvalues of Pn 
and Pn -|- Qn in {\z — A| < 5} coincide. 
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In other words, if we aim at detecting A with an error of it is 
enongh to consider estimations 

Pn{z) + Qn{z) =Z^- {2Mn - Fn)z + {[M^ - Gn) 

with snfficiently small ||-Fnil and ||Gn||, and hnd the spectrnm of Pn+Qn 
for snfficiently large n. The weights Wq and Wi are introdnced in order 
to allow independent control on how perturbations are measured in 
the truncations of M and For instance, two possibilities are: the 
absolute weights tco = wi = 1 and the relative weights 

^ AiWrW) ^ (/i/4+|A|)- 

3. Speed of convergence: examples 

In the examples presented below we hnd 5{n) explicitly for concrete 
operators M. Among other results, we illustrate how the error — A| 
depends strongly on the correct choice of We also provide numerical 
evidence suggesting that, in cases, the power 1/2 of 5{n) predicted by 
Theorem^ can actually be improved to a power a ~ 1. 

3.1. Finite rank perturbations of multiplication operators. As 

a hrst example we consider M = S + K acting on / G L^(—tt, tt) =: 
where 

i 

Sf{x) = s{x)f{x), Kf{x) = ^{f,gj)gj{x), 

i=i 

s(a;) is a bounded real-valued function and gj G L^. Since K has 
hnite rank, Weyl’s theorem ensures that SpeCess {P + K) = SpeCess S = 
Specs' = essRanges(a;). We aim at hnding asymptotic upper bounds 
for S{n) in the limit n —> oo in terms of invariants of s{x) and gj{x), for 
suitable £„ specihed later. The discrete Schrodinger operator studied 
in [H], [ini Example I] and the example discussed in (U Section 4], all 
satisfy the hypothesis of Lemma El below. 

Let Bi := By declaring the inner product normalized as 

ifyd) •= ^ clearly Bi is an orthonormal basis of L^. 

Denote by h{n) := (h, the Fourier coefficient of a complex-valued 
real function h. The following notation is in place for the examples of 
this section. For g > 1, 

0‘^ := {h : 27r — periodic in M s.t. \h{n) \ = 0(|n|“'^) as n —> ±cx)}. 

The well known rules relating Fourier coefficients and differentiation, 
easily show that if h' G then h G Hence, h G whenever 





NON-VARIATIONAL... 


7 


h G C"’(M) is 27r-periodic. Furthermore, h & for functions h of 
bounded variation, cf. e.g. U §2.3.6], 

Lemma 3. Let M = S + K and Cn = Span ..., If 

G for all j = 1,..., I, then we may choose 5{n) = 

as n ^ oo in 

Proof. Since M = S + K is bounded, it is enough estimating the left 
hand side of for k = 0. Indeed, 

WUnM^Unfj - - x’^f^w 

<i\\M^\\ + \X\^)\\U,,fj-fj\\. 

Now, since Mfj = Xip, then (A — s{x))ip{x) = ^ 9j)9j{^) and so 

the hypothesis guarantees that 'ip{x) G (9^. The estimate 

WUnip - V'lP = Y1 n ^ oo, 

|j|>n |j|>n 

completes the proof. 

In particular we may formulate the following conclusions for M = 
S + K and £„ = Span ..., e®®®^}. 

a) If s{x) and 9j{x) are 27r-periodic functions in C'®'(M), then we 
may choose 5{n) = 0{\n\~^~^/‘^) in ©• 

b) If s\x) and 9 'j{x) are of bounded variation, then we may choose 
5{n) = 0 (|n|“^/^) in ( 0 . 

c) If s(x) and 9j{x) are of bounded variation, then we may choose 
5{n) = 0 (|n|“^/^) in ( 0 ). 

All these assertions follow from the fact that if 
A G SpeCdisc (“S' + A'), then necessarily A ^ essRange s(x) and so 
is of the same degree of smoothness as s{x). The interesting case in 
our current discussion is c), because gaps in the essential spectrum only 
occur if s{x) is discontinuous. 


3.2. A numerical example. Let M = 5 + A be as in the previous 
section. Assume that 


f —2 + sin(2x), —TT < X < 0, 
I 2 + sin(2x), 0 < X < TT, 


I = 1 and 9 i{x) = y/2, so that Kf{x) = 2/(0). We may hnd SpecM 
in closed form. Clearly SpeCess Af = [-3,-1] U [1,3]. On the other 
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n 


log( z^-\) 

log(n) 

Slope 

190 

0.040879 

-3.1971 

5.247 

-0.50849 

235 

0.036691 

-3.3052 

5.4596 

-0.48708 

280 

0.033689 

-3.3906 

5.6348 

-0.50956 

325 

0.031226 

-3.4665 

5.7838 

-0.4876 

370 

0.029312 

-3.5297 

5.9135 

-0.50963 

415 

0.027647 

-3.5882 

6.0283 

-0.48835 

460 

0.026291 

-3.6385 

6.1312 

-0.50928 

505 

0.025071 

-3.6860 

6.2246 

-0.48918 

550 

0.024046 

-3.7278 

6.3099 

-0.50875 

595 

0.023103 

-3.7678 

6.3886 

-0.49003 

640 

0.022292 

-3.8035 

6.4615 

-0.50813 

685 

0.021535 

-3.8381 

6.5294 

-0.49086 

730 

0.020873 

-3.8693 

6.593 

-0.50748 

775 

0.020249 

-3.8997 

6.6529 

-0.49167 

820 

0.019695 

-3.9274 

6.7093 

-0.50682 

865 

0.019169 

-3.9545 

6.7627 

-0.49244 

910 

0.018696 

-3.9795 

6.8134 

-0.50617 

955 

1000 

0.018245 

0.017835 

-4.0039 

-4.0266 

6.8617 

6.9078 

-0.49318 


Table 1. 


Estimation of the exponent of 5{n) for ap¬ 


proximations of A_ using Bi- 


hand, according to ^ Lemma 7], we know that A is an eigenvalue in 
the discrete spectrum of M if and only if 


dx 


+ 


dx 


(A -I- 2) — sin(2x) Jq (A — 2) — sin(2x) 


= TT. 


The explicit computation of the integrals reveals that SpeCdisc M con¬ 


sists of two eigenvalues, the variational A_| 
prone A_ ~ —0.7674. 

A direct calculation yields 


s(j) = 


4 

ijn ’ 

(^ 2 , 


3.5796 and the pollution- 


j - odd. 




even, 


so that |s(j)| ~ as j — > ±cxd (here denotes the Kronecker 

symbol). Thus, by Lemma 0 we can choose 5{n) = and, 

according to Theorem [TJ the existence of G SpecP„ such that \z^ — 
-^±1 = in the limit n —>• cx) is predicted. 

From the explicit expression for ^(j) it is not difficult to find Pn{z). 
In Table ^ we report on the numerical approximations of \z~ — A_| 
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for different values of n = 190 : 1000 and the corresponding pairwise 
slopes between the steps n and n + 45 of the graph log \ z~ — \_\ vs 
log(n). We have found the numerical data by writing explicitly Mn 
and and computing the eigenvalue of the companion matrix of 

Pn{z) that is nearer to A_. For calculations we use the standard eigs 
routine available in the Matlab package. The last column strongly 
suggests that \z~ — A_| = 0(n““) for a ^ 1/2. Similar numbers are 
found for the variational A_|_. 

Table [T] suggests that there is generally a signihcant gap between 
the approximation predicted by & and the actual rate of convergence 
—A| —>■ 0. An obvious conjecture is that, perhaps, the bound 6{nYP 
in Theorem^ can actually be improved to S{n). 

3.3. Direct sum of multiplication operators. Choosing the right 
basis is crucial for achieving efficient approximation. In the example 
discussed above, this choice should be made attending the nature of the 
symbol, i.e. piecewise continuity. Below we introduce the correct basis 
to deal with symbols such as ©• The following results have obvious 
extensions to direct sum of any finite number of operators, i.e. many 
gaps in the essential spectrum. In order to keep our notation simple, 
we only consider two summing terms. 

Assume now that M = S' + iF is an operator acting on 
where 


( 10 ) 



Si, S 2 are real-valued bounded symbols and Gj = ((y'yi, G L^©L^. 
The difference with the previous model lies on the fact that SpeCg^g M = 
essRange si(a:) UessRange 52 ( 0 :), so gaps in the essential spectrum may 
occur even when si and S 2 are continuous and 27r-periodic. 

The following lemma is the natural adaptation of Lemma El to the 
present situation. Let e„(x) := and h„(x) = (0,e*"''^)L Then 

{cn, is an orthonormal basis of with the usual inner 

product. 

Lemma 4. Let M = S + K where S and K are as in (fTn|l . Let 
Cn := Span {e_„, h-n, • • • e^, h„}. If Y-frlil) ^ ^ = 1,..., / 

and m = 1,2, then we may choose 6{n) = as n ^ 00 in 

®. 

The proof is analogous to that of Lemma El 
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n 




log( - A 

) 

log(n) 

Slope 

12 


0.037578 


-3.2813 

2.4849 

-2.8383 

18 


0.011889 


-4.4322 

2.8904 

-3.1543 

24 

0.0047977 

-5.3396 

3.1781 

-4.9788 

30 

0.0015796 

-6.4506 

3.4012 

-4.9366 

36 

0.00064217 

-7.3507 

3.5835 

-7.194 

42 

0.00021185 

-8.4596 

3.7377 

-6.7381 

48 

8.6154e-05 

-9.3594 

3.8712 

-9.4147 

54 

2.8424e-05 

-10.468 

3.989 

-8.5396 

60 

1.156e-05 

-11.368 

4.0943 

-11.635 

66 

3.8138e-06 

-12.477 

4.1897 

-10.341 

72 

1.5509e-06 

-13.377 

4.2767 

-13.86 

78 

5.1145e-07 

-14.486 

4.3567 

-12.191 


Table 2. Estimation of A_ using B 2 


As an application of Lemma |31 we now consider a better suited basis 
for estimating the eigenvalues of the symbol discussed in Section Id.21 
The crucial observation is that is iso metrically isomorphic to 

via the unitary map 


fix) 




The orthonormal basis {e„, /i„} maps under U* onto an orthonormal 
basis B 2 = {lA*enM*hn} C in such manner that Wcn span all 
functions with support contained in [—vr, 0] and U*hn all those with 
support in [0,7r]. When s{x) is as in Q, 


(aw) 


/ (-2 + sina;)/i(a;)\ 
(2 + sina:)/ 2 (a:) ) 





Then, Lemma |3] along with Theorem ^predict that the basis B 2 would 
produce sequences approaching to A± whose speed of convergence is 
faster than any power of n. That is, if M is as in Section EUl and 


Cn = Span {U*e-n, U*h_n, • • •, Whn}, 

then for all g > 0 there is c{q) > 0 independent of n such that \z'f — 
A±| < c{q)n~'^ for all n G M. 

In Table 121 we report on computations of the nearest point in Spec Pn 
to A_ and list the values analogous to Table ^ In this case we only 
consider n = 12 : 78. Notice that the approximation of A_ using B 2 
at n = 12 is already more accurate than the one obtained using the 
basis i3i at n = 1000. The super-polynomial speed of approximation 
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predicted by Lemma 0] and Theorem ^ is already evidenced by the hrst 
few entries of the last column. 

3.4. Schrodinger operators with band gap essential spectrum. 

For our last example we consider an unbounded operator. Let M be 
the one-dimensional Schrodinger operator 

(11) Mf{x) = -dlf{x)+V{x)f{x), f G DomM := 

acting on L^(R), where the potential V = Vi + V 2 , Vi e and 

V 2 G Here denotes the derivative Sobolev 

space: V G if and only if H G Lp(R) and d^V G L^(M) for 

q= 

Since, in particular, V G L^(R) + L°°(R), multiplication by V is 
bounded with relative bound 0 [THl Theorem XIIL96]. Recall that the 
operator L is said to be Lq bounded, if and only if DomLo ^ DomL 
and 

\\Lf\\ < a\\Lof\\ + PWfW, f G DomLo 

for suitable non-negative constants a and jS. The inhmum of all a allow¬ 
ing the above for some (3 is called the relative bound of L (with respect 
to Lq). By virtue of Kato-Rellich theorem |17[ Theorem X.12], the 
above choice of DomM guarantees M = M*. Notice that DomM^ = 

Furthermore, since Vi is relatively compact with respect to —d^ 
Theorem XIIL15] and V 2 is bounded, Vi is also relatively compact 
with respect to —d^ -|- V 2 . Hence, by Weyl’s theorem, the essential 
spectrum of H is characterized completely by V2, that is SpeCgg^M = 
SpeCgss (—+ V 2 ). If 1^2 is periodic [IHl Chapter XHL16], the essential 
spectrum is bounded below, it lies in bands and it extends to inhnity. 
In general, Vi may produce non-empty discrete spectrum. 

We put 

(12) £„ = span{(^o, ■ ■ ■, 0n} 

where (f>j{x) = Ckhj{x)e~^^^‘^, cj = , , and hj{x) is the Hermite 

polynomial given by Rodriguez’s formula 

h,{x) = {-iye^\die-^^). 

The choice of Cj ensures that \\4>j\\ = 1 for all j = 0,1,.... It is 
well known that (pj{x) are the eigenvectors of the quantum mechanical 
harmonic oscillator H := —-|- so that is an orthonormal 

basis of L^(M). 
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Theorem 5. Let M be given by (HH) and Cn be given by m- Then for 
all X G SpeCjjg^M there is Zn G Spec Pn such that Zn —>■ X as n —>■ oo. 
Furthermore, if d^V (t) is continuous and bounded for some q > 3, 
then \zn — X\ = o{n^) as n ^ oo. 

The proof of this result will be a consequence of three technical 
lemmas. 

Lemma 6. Ifd^V{x) is continuous and bounded for some q>l, then 
tf G and x G M, for suitable constants 

di,ai > 0. 


Proof. Since V G then V'(x) is a continuous function. By 

hypothesis, if" = [V — X)if, showing that if"{x) is also continuous. By 
repeating recursively this argument q times, the equality 

(13) = i-diiV -\) + ...+ {V- \)dii, 

implies the existence and the needed continuity of 

By virtue of the results of m §C.3], the hypothesis on V ensures 
that 

(14) |i/i(x)| < xgM, 

for some positive constants ^ 2 , 02 - Since V{x) — X is continuous and 
bounded, then \if"{x)\ < 1 / 26 ““^'^'. By repeating recursively this argu¬ 
ment q times, (fTtlll and the hypothesis, yield the desired upper bound 
for 


Lemma 7. Let M be given by 03 . Then any eigenfunction of M lies 
in Dom . 

Proof. From the inclusion if G DomM^ = hF^’^(M) and (fTT|l (which 
do not require any smoothness property on V), it follows that if G 
hF^’^(M)nDom (x^) = Domif and {—if"+x‘^if) G DomH. This ensures 
the desired property. The only non-trivial facts in the latter assertions 
are, perhaps, the inclusions x'^if G fF^’^(M) and if" G Dom(x^). The 
hrst follows from the second, by differentiating twice the term x'^if and 
noticing that {xif')' G L^(M). The second one is achieved as follows. 
Since Mif is an eigenvector of M, | — if"{x) -|- V{x)if{x)\ < 
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Then 


< (i4 + 


\ 1/2 

x'^\V{x)\^ \'ip{x)\^(\x ) 




1/4 T c/o 


x4e-2“2N|l^(x)|2dx) 


\ 1/2 


The latter integral is bounded because of Id G L^(M) + L“(R), therefore 
-ifj” G Dom (x^). 


Lemma 8. Let M be given by (HID. Then both M and are 
bounded. 

Proof. Multiplication by Id is d"^ bounded with bound 0. Then M is 
bounded with relative bound 0. Similarly, from the identity 

{-dl + Vfu = d^u - 2Vdlu - 2V'd^u + {V^ - V")u 

and the fact that V, V', V" lie in L^(M) + L“(M), one may deduce that 
is 94 bounded with relative bound 1. Then the proof reduces to 
showing that is bounded. For this, let 

A-.= 2-^/‘^{x + d^) and A* = 2-^/‘^{x - d^). 

Then (—9^ + x^) = 2(74^4* — 1) and = {A — A*Y. Thus the desired 
property follows from the identity fcf. [T71 eq.(X,28)|) 

II^# 1 ... A*’^u\\ < c||(—+ x^)'^'^^m||, g = 1,2,... 

where A^'"^ is either A or A*. This is easily shown by induction and 
using the estimate 

\\{-dl + xrfM < wi-dl + x^yi^w, p<q. 

This completes the proof of the lemma. 

We may now complete the proof of Theorem HJ Let 0 := The 

crucial point lies in the following observation. According to Lemma |H1 
there exist non-negative constants a and /3, such that 

- A"0|| < - a^II 

<||M'^(n„0-0)11+ ||A'^(n,,0-0)11 
< «i|Lf^(n„0 - 0)11 + (0 + |A|^)||n„0 - 011 
= a||n„0-0|| + (0 + |A|'^)||n„0-0||. 

This ensures directly the hrst part of the theorem. For the second part 
we only require estimating ||n„0 — 0|| and ||n„0 — 0|| as n —> cx). 
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A straightforward application of lemma IHl shows that, if 
d^V{x) is continuous and bounded for q > 2, then for any set of 
polynomials {po{x),... ,Pq+ 2 {x)}, the functions Yl‘j=oPj{^)4’^^K^) 
Yl'jZo ci-re continuous and square integrable. By virtue of 

the fundamental relation 


hj (x) 


h'j+i{x) 

2(J + 1)’ 






integration by parts yields 


{'tp,(t)k) = Ck I 'ilj{x)hk{x)e = Cfc / i^{x) 

J — oo 

[^p{x)e~'"^/‘^\hk+i{x) dx 


' —CO 


\ -.V2 d. 


2{k + 1 ) 

'2{k+ !)]_, 


CO 

OO 


[i)'{x) — x'ip{x)]e ^^^"^hk+iix) dx = ... 


-ly^'^Ck 


2'?+2(fc + 1) ■ ■ ■ (fc + g + 2) 

(-1)'?+^(C, 0fc+g+2) 

2^ ^J{k+ 1) ■ ■ ■ {k + q + 2) 


q -\-2 


'^Pj{x)'ilj^^\x) e '^^/‘^hk+q+2{.x)dx 

3=0 


where is continuous and square integrable. 

Thus 


k=n-\-l 

^ _ 1(C) 0fc+g+2)P _ 

- ^2‘'+2(fc + l)---(A: + g + 2) 

^ 2,+2^^g+2 l(C,0A)+,+2)P = o(n-(''+2)) 

as n —>■ cx). A similar calculation with 0 instead of 0, gives 
||n „0 - 011 ^ = n ^ CX), 

completing the proof of Theorem El 


We remark that Theorem El may be extended to higher dimension 
without much effort. 
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4. Distance to singularity, the spectral function 


In order to simplify the notation in the proof of Theorems Q and 
121 we hrst discusses the notion of distance from P{z) to the nearest 
singular matrix. For convenience, all the ideas in this section refer 
to general matrix polynomials of arbitrary degree P{z) = 
where Ak G 0^^ and detA^ ^ 0. Some of the results below are 
related to the recent electronic manuscript by Davies [7j and the paper 
by Lancaster and Psarrakos |14j . 

For all z G C, we dehne 


ap{z) 


inf 

v^O 


\\P(z)v 



Four different characterizations of this function are, 

a) ap{z) equals the smallest singular value of P{z), 

b) ap{z)~^ = ||P( 2 ;)“^|| whenever det P( 2 ;) ^ 0, 

c) ap{z)-^ = sup„„^o whenever detP(^) ^ 0, 

. apiz) = min{||P|| : det[P(; 2 ) + P] = 0, P G 

^ = min{||P|| : det[P( 2 ;) + P] = 0, RankP = 1}. 

Since 

SpecP = {z G C : ap{z) = 0}, 

we may call the scalar quantity cp^z) a “spectral function” of the 
matrix polynomial P. 

The proofs of a), b) and c) are straightforward and property d) 
holds trivially when ap{z) = 0. In order to prove d) when ap{z) ^ 0, 
we use b). If ||P|| < ||P( 2 ;)“^||“^, then ||P( 2 ;)“^P|| < 1 and P{z) + 
P = P{z){I + P{z)~^E) so that det(P( 2 :) + P) 7 ^ 0. This shows that 
||P(z)“^||~^ can not be greater than the right sides of d). Conversely, let 
V G 0 be such that ||P(z)“^n|| = ||P(z)“^|| and ||n|| = 1, and let u : = 
\\P{z)~^\\~^P{z)~^v. Then ||m|| = 1. Put Ew := —||P( 2 :)“^||“^(tc,M)n 
for all w G CL Then the linear operator P has rank one, it satishes 
||F;|| = ||P(z)-i||-i and Eu = -||P(z)-i||-i'i;. Thus (P(z) + P)m = 0 
so that det[P( 2 ;) + P] =0. This ensures that the bottom right side of 

d) is not greater than ||P( 2 ;)“^||“^. 


The following lemma will play a fundamental role in the sequel. See 
in Lemma 2] and [71 Section 3.4]. 


Lemma 9. The non-negative function ap{z) is Lipschitz continuous 
in every compact subset of the complex plane. Furthermore ap{z)~^ is 
subharmonic in C\ Spec P, so that zq is a local minimum of cp if and 
only if zq G SpecP. 
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Proof. The first part of the lemma follows easily from the triangle in¬ 
equality once we have proven the estimate 

(15) \ap{z) — ap{w)\ < \\P{z) — P{w)\\ for all z, w G C. 

We show the latter by considering two separated cases. If 2 ; G Spec P 
and w ^ SpecP, since P{z) = P{w) + {P{z) — P{w)) is not invertible, 
d) ensures If z,w ^ SpecP, then 

P{z)-^ - P{w)-^ = P{z)-^[P{w) - P{z)]P{w)-\ 

By the triangle inequality, 

||PW-‘|| - ||PW-'||| < ||PW-'||||P(tt.) - PW||||P(tt,)-'|| 

and so m is consequence of b). The second part of the lemma follows 
from c) and the elementary properties of subharmonic functions. 

We now consider weighted perturbations of P{z) in the sense studied 
recently in uni and M- We will require these results in the proof of 
Theorem 121 Below, Q(z) = shall always refer to a small 

perturbation of P. Let £ > 0 and w = (tco, • • •, Wn) where Wk > 0. The 
(weighted) e-pseudospectra of P are the set given by 

_ p f „ p (p . det[P(z) -L Q{z)] = 0, for some Q{z) = EkZ^ f 
' I ' such that llPfcII < ^ = 0,..., m /’ 

This dehnition was studied by Higham and Tisseur in PI and it 
extends the standard dehnition of pseudospectra of a matrix A in the 
obvious manner, by considering P{z) = {z — A) and w = (1,0). The 
weight w is introduced in order to allow freedom in controlling the 
perturbation of each individual coefficient of P, e.g. this may be given 
in an absolute sense {wk = 1 ) or in a relative sense {wk = ||^fc||)- 
Observe that Speco,^P = Spec P. 

By combining the remarkable result m Lemma 2.1] with b), we 
achieve the useful characterization 

Spec£,ujP = {z e C : ap(z) < e(wo + wi\z\ + ... -1- Wm\z\'^)}. 

Notice that we do not require ||Pm|| to be small, so Spec^^wP is not 
guaranteed to be bounded. This shall not be an important point here. 
Necessary and sufficient conditions for these set to be bounded may be 
found in [TH Theorem 2.2]. 

The following lemma is relevant in the proof of Theorem |21 This has 
also been established by Lancaster and Psarrakos in the more sophis¬ 
ticated m Theorem 2.3]. 

Lemma 10. Let e > 0. Let PL he a connected component of SpeCe^u, P 
such that n n SpecP 7 ^ 0 . If ||Pfc|| < WkS for all k = 1 ,... ,m, then 
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f2 n Spec {P + Q) 7 ^ 0. Furthermore P + Q has the same number of 
eigenvalues in counting multiplicity, than P has. 

Proof. Let Ps{z) := + hEk)z^ for 0 < 5 < 1. Then Pq{z) = 

P{z) and Pi{z) = P{z) +Q{z). Since these eigenvalues are the zeros of 
a certain family of scalar polynomials whose coefficients depend con¬ 
tinuously upon S, the eigenvalues of Ps{z) depend continuously on S. 
Let fas G SpecPs, be such that po G hi. Since ||(5L^fc|| < \\Ek\\ < WkS, 
then fis G Spec^^^P for all 0 < 5 < 1. Being a continuous trajectory, 
necessarily G hi for all 0 < 5 < 1. 

5. Proof of theorems [U and [2] 

Throughout this section we shall write = ap„. 

Let 6i:=2-|-2|A|-|-3|Ap. By virtue of the triangle inequality, 

||n„(A - M)2n„v^|| = - 2XUnMUnf: + 

< WUnM^UniJ - X^UnfjW + px'^Unfj - 2XUnMUnij\\ 

< ||n„M2n„V' - AVil + 2|A| IIAV^ - n„Mn„i^|| + 3 |aH|v^ - 

< {b,-l/2)6{n). 

Then, from the definition of cr„, by choosing v = Hn'f, we achieve the 
following. 

Lemma 11. There exists Ni > 0 such that a„(A) < bi6{n) for all 
n > Ni. 


The proof of both theorems depends heavily on applying the second 
part of Lemma El For this we combine Lemma im with lower bounds 
for an{z) along the boundary of small neighbourhoods of A. Our first 
step towards computing these bounds, involves finding estimates for 
||-Pn(^)~^|| in subspaces that are orthogonal to the eigenspace associated 
to A. 

Let £ := span {f) G DomM : Mfj = X'lf}, d\ := dimT and Ll^- be the 
orthogonal projection onto £. Let /i = ±p be such that A —/i G Spec M 
and 

A := M(J — Ilf) 4- (A — fi)T\.£, DomA = DomM. 

Then, by construction, A = A* and Spec A = SpecM\ {A}. For z G C, 
we write A{z) := {z — A^ with Dom74(z) = DomM^, K{z) = fi{f2\ — 
2z — so that 

A{z)x = {z — MYx — K{z)x for x G DomM^, 

An{z) := I\.nA{z)\Cri, Kn{z) = Il„K{z)\Cn and An{z)-^ = [An{z)]-'^. 
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Let D ■= {\z — \\ < jj,/A}. Then D does not intersect Spec A and A 
is the closest point in SpecM to the boundary of D. Furthermore, if 
z E D, then 

inf Ref^: —inf (w — A — Re^:)^ — 

^,A+^] it;GR\[A—^,A+//] 

> /iV2 > 0. 

Put 7 := /i^/2. Since — A)^|£„n,n) = {{z — A^v^v) for all 

V E Cn and A{z) is a normal operator, then the following inclusion for 
the numerical range of An{z) holds for all G -D, 

NumR„(;2) C NumR(2;) C Conv[Spec {z — Ay] 

C {[z — wY : w E [A — p, A + n]} C {Re {z) > 'j > 0}. 

Thus An{z) is invertible and 

||Al;^(^)|| < [dist (0,NumAl„(^))]“^ < 7“^ 

for all 2: G -D and n E N. 

Let 

62 •= 1 + 2(/i/4 + |A|) + (/i/4 + |A|)^ + 2/i(4|A| + 3/i/2). 

The triangle inequality ensures that 
\\An{z)Ilny - A{z)y\\ 

= ||n„(2: - - (z - Mfy - /i(2A -2z- ^j)[y - n„n£-n„i/]|| 

< - \^y\\ + 2\z\\\IinMIiny - Ai/I| + 

(2|/i(2A -2z- /i)| + |^P)||i/’ - n„i/|| 

< (1 + 2\z\ + {zl"^ + 2/i|2A — 2z — /i|)(5(n) < b 2 S{n) 

for all 2; G -D and n E N. Let u := {\ — fi — zY, so that A^z)^ = 

Then 

\\An{z)~^Tlny - A{z)~^y\\ = - z/“Vll 

< \\An{z)~^Iini) - + ||z/“^n„i/ - Z/“V|I 

< - An{z)Y\.ny\\ + |z/|“^||n„i/ - y\\ 

= 7"Vrln„A(^)i/ - An{z)iiny\\ + |z/|“^||n„i/ - y\\ 

< 7“Vrl^(^)^ - An{z)Uny\\ + |z/|“^||n„i/ - y\\ 

< (7"^&2 + 

< (4(7"^&2 + l)/i/3)5(n) =; h6{n) 
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for aW z ^ D and n E N. Hence, a straightforward computation shows 
that 

(16) \\An{z)-^UnK{z) - A{z)-^K{z)\\ < dxh5{n) 

for aX\ z E D and n G N. 

The following estimate, which is a direct consequence of the dehni- 
tion of K{z), will be needed below: 

||-^(^)|| < h(2|:2 — A| + /i) < 3/i^/2 =: 64 for all z E D. 


5.1. Proof of Theorem 12 By virtue of LemmalHl the only local min¬ 
ima of (Jn{z) are those points in SpecP„. Then the proof of Theorem^] 
reduces to hnding 6 > 0 and N 2 > 0 , both independent of 2 ; and n, 
such that 
(17) 

cT„(z) > cr„(A) for all 2: E {\z — A| = and n > N2. 

Since Spec Pn is hnite, the existence of 6 > 0 as we require in the thesis 
part of the theorem becomes obvious. 

Inequality (jl7|l is a consequence of Lemma and the following. 

Lemma 12. let b := (4&4 max o,nd let N2 E N be 

such that 


6 {n)<mm{b ^(/x/ 4 )^ (IdA^s) n > N2. 


Then crn(z) > biS(n) for all |;2 — A| = and n > N2. 


Proof. Throughout the proof we assume that n > N2 and I 2 ; — A| = 
Since < /i/4, then A{z) -|- K{z) = {z — M)^ is 

invertible and 


(18) 


\\[I+A{z)-^K{z)]-^ = \\{z-M)-^A{z)\\ 

= \\{z-M)-^[{z-Mf-K{z)]\\ 

= \\I -{z- M)-^K{z)\\ < 1 + b 4 {z - M)-^\\ 

< 1 + 64 [clist {z, Spec M)]-^ = ^ \ =■ [c(n)]“^ 

b^d[n) 


Thus, given v E Cn, 
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Since 6^ > Ab^dxb^ and 6 {n) < {Adxb^) \ then 6^ > iJ^s(n)^dxb3 

(uni) and an easy calculation yield ||An(^) ^^nK{z) — A[z) ^K{z)\\ < 

c(n)/2. Hence 

^^i|u|| < \\v + An{z)~^IinK{z)v\\ 

< \\An{z)~^\\\\An{z)v + IlnK{z)v\\ 

<r^\\A^{z-Mf\Cnv\\. 

Therefore An{z — is invertible and 

ajz) = ||(n„(s - M)2|£„)-I|ri > AA, 

Since 6^ > 7“^46461 and 5 {n) < 7(461)“^, it is easy to see that 6^ > 
■y-2biS(n ) ' ^lius cr„(z) > biS{n) as claimed in the lemma. 

This completes the proof of Theorem ^ Explicit expressions for b in 
terms of A and fi might be useful in applications. A direct substitution 
yields b = (max {65,126i} + 1)^/^ where 

, , 16 + 16|A|2 + 8/i + 57/i2 + 8|A|(4 + 17/i) 

05 dx -. 

/i 

5.2. Proof of Theorem m The proof will be a consequence of lem¬ 
mas MM and the following. 

Lemma 13. For all 0 < 6 < /i/4, there exists > t) independent of 
z or n, sueh that 

(5^7 

> WTA 

for all n > N 3 and <5 < |^ — A| < /i/4. 

Proof. Let C := {5 < 1^; — A| < /i/4}. Throughout the proof we assume 
z e C. Since C C D, then A{z) + K{z) = {z — M)^ is invertible. By 
virtue of the third line of m, 

II [/ -h A{z)~^K{z)]~^\\ < 1 64[dist {z, SpecM)]”^ < ^ =■ cf^. 

Let A"3 > 0 be such that dxb^d^n) < C2/2 for all n > N^. Then a 
straightforward computation along with (uni) yield 

\\A„(z)-^n„K(z) - A(z)-'K(z)\\ < ^ 

for all n > N 3 . Given v E Cn, then 
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Figure 1. Proof of Theorem |21 The shaded region rep¬ 
resents Spec^^Pn which is symmetric with respect to 
the real axis (dotted line). The dark intervals represent 
SpeCgggM and the cross the eigenvalue A. 

Thus C 2 i|ni |/2 < \\Pn{z)v\\'j~^, so that Pniz) is invertible and (Tn{z) > 
7 C 2/2 for all n > N 3 . 

We now complete the proof of Theorem |21 Notice that the case £ = 0 
is an obvious corollary of Theorem ^ Let £ > 0 be as in the hypothesis 
of Theorem 121 and let w = (tCojiTijO). A straightforward argument 
yields 

D n SpeCe^^s Pn = {z e D : (T„(z) < e{wo + Wi|^|)} 

C {. 6 C : 

Thus, if n > N 3 , by virtue of Lemma UHl C fl Spec^^u, Pn = 0 guaran¬ 
teeing the hrst conclusion of the claimed result. 

On the other hand, Lemma im ensures the existence of A 4 > Ni 

c2 

such that (T„(A) < 2 { 5 ' 2 +b 4 ) ^ ^ iLiax{A^ 3 , A 4 } and 

assume that n > N. By Lemma El there exists G Spec Pn such that 
— A| <6. Furthermore, the connected component fin C SpeCs^w Pn 
such that Zn & fin, satishes fin C D\C. Hence, the second and third 
conclusions of the theorem follow directly from Lemma E3 

6. The essential spectrum 

The essential spectrum of M is usually found by means of analytical 
methods. Nonetheless, besides of being a natural question per se, the 
numerical evidence in all the examples discussed in Q. US. i and 
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[S], suggest that approximation of this portion of the spectrum may 
also occur in the “second-order” projection method described above. 
In this hnal section, we discuss some results and open questions related 
to this issue. 

Ideally, we would like to know where does the whole set Spec Pn 
accumulates in the limit n —> oo. To this end, we may consider the 
following two limiting set and then study the connection between them 
as well as their relationship with Spec M. Given e: > 0, let 


Ag := {C G C : there exists Zn ^ ( such that cr„(^„) < e} 
and := {C G C : lirnsup^^^ cr„(C) < e}. 

Theorem Q ensures that Specdisc M C M fi Aq. We now ask, which 
conditions yield a similar property for SpeCess M. We do not intend 
to answer this question here. Nevertheless, elementary properties of 
these set might provide an insight towards further investigations in this 
direction. Notice that A^ is the (uniform) limit set of SpeCe,(i,o,o)-Pn as 
n —> cx). 

For simplicity we assume from now on that M is bounded. Thus 
(19) Aq C So = [^ A^. 

5>0 5>0 


Indeed, directly from the dehnition, it follows that 

\Pn{w)v\\ + 1^ — tc|||n„(2M — z — t(;)n„u| 


an{z) < min ■ 


( 20 ) 


< <Jn{w) + \w — z\ sup 

< an{w) -f |tc — z\c 


II {2Mn — z — w)v\ 


for all |tc — z| < £ < 1, where c > 0 is chosen independently from e, n, z 
and w, because of M is bounded. Now, assume on the one hand that 
2 ; ^ Sg. Then there exists n{j) G N and a > 0, such that an(j)(z) > e+a 
for all j G N. According to (HD, 


an(j){w) > an(j)(z) — c\w — z\ > e P a — ce 

for all \w — z\ < e. Hence, by choosing e = a/{2c), it becomes evident 
that an(j){w) > e whenever \w — z\ < e for all j G M, so that 2 : ^ A^. 


Thus 


A. C S, 


for all £ > 0. 


On the other hand, it is not difficult to prove that 


Sg C Ag +5 for all £ > 0 and 5 > 0 
and that S^ C Sq. These three inclusions ensure flTTHl . 




NON-VARIATIONAL... 


23 


The following proposition is crucial in the method suggested recently 
by Davies and Plum in [S]. 

Proposition 14. If M is bounded, then Sq fl M = SpecM. 


Proof. Indeed, if A G Spec M, for each /c > 0 there is fjk G H, W'lpkW = 1, 
such that II(A — M)‘^‘ipk\\ < I/k. Then 

lim < i/k 

n^oo ||n„l/>fc|| 

and so cr„(A) —>• 0 as n —> cx). Conversely, notice that if A G M but 
A ^ SpecM, then NumP„(A) C Num(A — M)^ C [fi,oo). Hence 

c’'n(A) > dist [0,NumP„(A)] > fi, 

so that A ^ So n M. 


In other words, the inclusion Sq C Aq complementary to the hrst one 
in (fT^. will automatically imply approximation to the whole spectrum, 
in particular to SpeCess M. For instance, Proposition 3], if M^ = Id, 
then So = Aq C {|C| = 1 }. The validity of this inclusion is closely 
related to the problem of whether there is an upper bound independent 
of n for the size of the blocks in the Jordan canonical form of An, 
in the factorization Pn{z) = [z — An){z — A*). Indeed, let Rn{z) = 
{z - Sn){z - S'*), where 


1 0 


\ 


G C 


nxn 


V 1 0/ 

Then SpecP„ = SpecS'n U Spec S'* = {0} for all n G N. Thus Ao = 
{0}. By choosing n = {z"'~^,... ,z,iy, Rn{z)v = 0,..., 0) so 

(^Rniz) < Izl”’*'^. Hence Sq = {| 2 :| < 1}. Of course, although strict 
inclusion holds in this case, it not clear whether Rn{z) = Pn{z) for 
some M = M* and {n„}. 


Acknowledgement. The author wishes to express his gratitude to 
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